AD-A203  480 


nun  Hit  '•'h 


TECH.  MKMO 
AMRO  2117 


UNLIMITED 


TCP],  . 

* '  *- 


tfjj  fa' ri$ 
y/nkjrv'W 

^iUKv 

4-14 

ROYAL  AEROSPACE  ESTABLISHMENT 


NUMERICAL  EVALUATION  OF  A  SOLUTION  REPRESENTING 
A  JET  IN  CROSS-FLOW 


C.  C.  Lytton 
J i  H.  B .  Smith 


August  1988 


1  i  LSj  ■  - ; 


DTIC 

TyiLECTE 

>Vl 

V  «  f  f  ft  5n?Q 


I 


Procurement  Executive,  Ministry  of  Defence 
Farnborouflh.  Hants 


UM.IM1VVD 


89  2  i  03: 


f 


i 

i 


i 


ROYAL  AEROSPACE  ESTABLISHMENT 
Technical  Memorandum  Aero  2137 
Received  for  printing  4  August  1988 

NUMERICAL  EVALUATION  OF  A  SOLUTION  REPRESENTING 
A  JET  IN  CROSS-FLOW 

by 

C.  C.  Lytton 
J.  H.  B.  Smith 


!  > 


SUMMARY 

A  recent  study  by  Needham,  Riley  and  Smith  of  a  jet  emerging  from  a  circular 
pipe  into  a  weak  cross-flow  under  the  assumptions  of  an  inviscid  and  incompress¬ 
ible  fluid  led  to  expressions  for  the  velocity  potential  in  different  regions 
near  the  orifice  as  definite  integrals  and  infinite  sums.  In  the  present  paper 
the  expressions  relating  to  flow  inside  the  pipe  and  jet  are  evaluated  with  great 
care.  It  is  confirmed  that  the  solutions  in  the  pipe  and  jet  are  continuous  across 
the  orifice  plane.  The  behaviour  of  the  solution  near  the  orifice  is  displayed 
in  graphical  form,  with  emphasis  on  the  consequences  of  the  failure  to  satisfy 
a  Kutta  condition  of  smooth  outflow  at  the  lip. 
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1  INTRODUCTION 

Needham,  Riley  and  Smith1  present  an  analysis  of  an  idealized  form  of  the 
jet-in-crossflow  problem.  The  flow  is  assumed  to  be  inviscid  and  incompressible. 
A  jet  emerges  from  a  semi-infinite  pipe  in  the  shape  of  a  circular  cylinder  into 
an  ambient  flow  whose  component  normal  to  the  pipe  axis  is  small  compared  with 
the  jet  speed.  The  configuration  is  illustrated  in  Fig  1.  The  flow  is 
irrotational  both  inside  and  outside  the  jet  boundary,  which  forms  a  vortex 
sheet.  This  sheet  is  assumed  to  lie  close  to  the  circular  cylinder  defined  by 
extending  the  pipe,  so  that  the  conditions  of  tangential  flow  and  continuity  of 
pressure  which  are  satisfied  on  the  jet  boundary  can  be  transferred  to  the 
cylinder . 

The  solution  obtained  is  in  two  parts,  corresponding  to  different  length 
scales.  If  the  ratio  of  the  cross-flow  velocity  component  to  the  jet  speed  is 
e  «  1  ,  and  the  pipe  radius  is  a  ,  the  first  part  of  the  solution,  correspond¬ 
ing  to  a  length  scale  of  a/e  ,  is  obtained  as  the  first  few  terms  of  a  power 
series  in  ex  ,  where  x  »  x/a  and  x  is  the  distance  from  the  exit  plane. 

The  second  part,  on  a  length  scale  of  a  ,  matches  the  first  part  to  conditions 
at  the  pipe  exit.  This  is  obtained  by  the  Wiener-Hopf  technique,  which  leads 
to  involved  expressions  for  the  velocity  potential  in  distinct  regions,  vis 
inside  the  pipe,  inside  the  jet,  outside  the  jet  and  outside  the  pipe. 

These  expressions  are  so  elaborate  that  there  is  no  obvious  way  to  confirm 
analytically  that  they  satisfy  the  required  conditions  of  continuity  across  the 
boundaries  of  the  regions.  Nor  is  it  possible  to  deduce  any  quantitative 
information  about  the  solution  by  inspection.  Since  there  is  at  present  very 
little  known,  in  a  mathematical  sense,  about  the  behaviour  of  jets  in  cross-flow, 
it  seems  worth  investigating  the  behaviour  of  this  model,  physically  unrealistic 
though  it  be.  Accordingly,  an  attempt  has  been  made  to  evaluate  numerically  the 
expressions  obtained,  first  to  confirm  that  the  appropriate  quantities  are  indeed 
continuous,  and  then  to  display  the  quantitative  nature  of  the  solution  near  to 
the  pipe  orifice. 

It  emerges  that  the  numerical  evaluation  is  far  from  straightforward;  the 
typical  task  is  to  evaluate  an  integral,  over  an  infinite  range,  of  products  of 
quasi-per iodic  functions  (changing  sign  infinitely  often)  and  a  function  which 
is  itself  defined  by  a  further  integration.  For  the  particular  case  in  which  the 
coflowing  component  of  the  ambient  stream  is  equal  to  the  jet  velocity,  ie  X  ■  1 
in  Fig  1,  one  of  the  quasi-periodic  functions  reduces  to  a  form  which  is  more 
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amenable  co  numerical  integration,  and  so  the  present  exercise  has  been  confined 
to  this  case.  The  restriction  to  X  »  1  does  remove  some  of  the  physical  con¬ 
tent  of  the  problem,  in  particular  implying  that  the  inclination  of  the  jet  to 
the  ambient  stream,  tan  *  (e/X)  ,  is  small.  Nonetheless,  the  particular  case  is 
believed  to  provide  an  adequate  test  of  the  analysis  and  a  useful  indication  of 
the  behaviour  of  the  solution.  Even  for  this  case,  considerable  care  and 
appreciable  computing  time  are  required,  so  it  is  thought  worthwhile  to  place 
the  methods  and  the  results  on  record. 

2  SUMMARY  OF  THE  ANALYTIC  SOLUTION  NEAR  THE  ORIFICE 

In  Ref  1  is  discussed  the  situation  of  an  incompressible,  inviscid  fluid 
emerging  with  unit  average  speed  from  a  long  straight  pipe  of  circular  cross- 
section  with  unit  radius.  (So  velocities  and  distances  are  scaled  on  pipe  flow 
speed  and  pipe  radius.)  The  jet  emerges  into  a  stream  of  similar  fluid  which 
has  undisturbed  velocity  components  X  parallel  to  the  pipe  and  e  perpendicular 
to  it,  as  indicated  in  Fig  t.  Cylindrical  polar  coordinates  (r,e,x)  are  taken 
aligned  with  the  pipe,  with  the  origin  at  the  centre  of  the  pipe  orifice, 
x  directed  out  of  the  pipe  and  0  measured  from  the  pipe  leeward  generator. 

Thus  the  pipe  surface  is  given  by  r  •  1  ,  x  {  0  . 

Since  the  cross-flow  parameter  e  is  small,  it  is  assumed  that  near  the 
orifice  the  extension  of  the  pipe  (r  «  1,  x  >,  0)  can  be  used  to  separate  the 
representations  of  the  flow.  In  equation  (3-23)  of  Ref  1  the  velocity  potential 
$  in  r  >  1  is  expanded  formally  in  powers  of  e  .  This  expansion  may  be 
written  as 

<j>  »  Xx  +  e^r  +  -^cos  9  +  tZ  ^gCr.x)  +  i)j2(r,x)cos  26 J  +  O^e3)  ,  (1) 

and  the  potential  $  in  r  .$  1  is  similarly  expanded  as: 

$  ■  x  +  e2^i|i0(r,x)  +  i(/2(r,x)cos  20 

Here  ^  ,  i|/q  are  functions  giving  the  axially  synmetric  part  of  the  solution, 

which  will  not  concern  us  further,  and  ^  8ive  the  part  of  the  solution 

depending  on  9  ,  to  this  second  order.  In  a  small  change  of  notation  from  Ref  1, 

a  factor  cos  20  has  been  taken  out. 
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For  the  region  inside  the  pipe,  x  4  0  and  r  ^  1  ,  the  solution  for  # 
is  obtained  from  equation  (3-36)  of  Ref  1  as: 


72(r«x) 


/  , 

V  1  +  X  /  K+(ifln)J2 W 


where  the  sunmation  is  over  the  positive  zeros  to  of  dJ^Cpi/dp  S  j^(p)  “  0  , 
where  J2  is  the  Bessel  function  of  the  first  kind  and  second  order.  To  define 
K+  ,  we  first  introduce  a  function  y(s)  of  the  complex  variable  s  : 


(s  ♦  «)*<•  -  «>* 


where  branch  cuts  are  taken  from  ±id  to  ±i°»  ,  and  the  square-root  branches  are 
chosen  to  have  positive  real  parts  for  s  real.  Hence,  when  s  «  u  is  real, 

Y  “  | in |  .  Then,  from  the  equations  (B-13),  (B-16)  and  (B-15)  of  Ref  1: 


VY)  .2 


k2(y)  I 

KpTTj 


F±(s)  -  ± 


’  [  — — —  Urn 

2,1  L  )  +  A2  J 


('  *  »’)  iS 


(vs)! 


where  the  upper  or  lower  signs  are  to  be  taken  throughout.  K_  is  needed  later. 
Here  I2  and  are  modified  Bessel  functions  of  the  first  and  second  kinds 

respectively,  and  C+  ,  C_  are  contours  of  integration  lying  respectively 
below  and  above  all  singularities  in  the  complex  u-plane  .  In  the  following 
section  it  emerges  that  the  contours  can  be  taken  to  lie  along  the  real  axis. 
Finally,  A  is  given  as  the  quotient  of  two  integrals: 
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where 


l  &  .iK(p).J'(p).K+(ip)  , 
0 

iK(p) . J’(p) .K+(ip)  + 


(8) 


(9) 


in  which,  since  p  is  real  and  positive,  we  have  from  (B-19)  of  Ref  1 : 


iK(p) 


iK^(-ip) 

D(-ip) 


(10) 


where  Re  denotes  the  real  part  of  a  complex  function,  and  for  convenience,  we 
have  introduced: 


D(-ip)  -  I2(-ip)K’(-ip>  -  X2K2(-ip)l’(-ip)  .  (11) 

We  have  used  I2(ip)  *  I2(-ip)  and  I2(ip)  “  -  I2(-ip) ,  see  Ref  2. 

For  the  region  inside  the  jet,  :  )  0  and  r  4  1  ,  the  solution  for 
is  obtained  from  equation  (3-38)  of  Ref  1,  after  correcting  a  misprint  in  the 
sign  of  Ap  ,  as: 


^(r.x) 


1  +  X 


tt/2  ^1  +  X2)  0 


CO 

f  f*E 

/  2 


iK(p) . J2  (pr)  .K_(-ip).e 


-px 


.(1  -  Ap) 


(12) 


For  the  region  outside  the  jet,  x  )0  and  r  1  ,  the  solution  for  is 

obtained  from  equation  (3-42)  of  Ref  1  as: 


xp2<r,x) 


Xx 


0  +  w/2(l  ♦  X2)* 


.2  Re 


K2(-ipr) 

D(-ip) 


•  K_("ip)  .  iljdp)  .e  px.(1  -  Ap) 


(13) 
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For  the  region  outside  the  pipe,  x  ■(  0  and  r  >,  1  ,  the  solution  for  is 

obtained  from  equation  (3-40)  of  Ref  t  as: 


Cr- 


s/I^I  +  X2) 


.2  Re 


K2(-ipr) 

K'(-ip) 


e~p|*l 

K+(ip) 


. (1  +  Ap) 


04) 


3  PRELIMINARY  MANIPULATIONS 

The  Bessel  functions,  and  modified  Bessel  functions,  of  orders  0  and  1 
with  real  and  positive  arguments,  are  available  as  subroutines  in  the  NAG  Fortran 
Library.  To  obtain  the  functions  of  order  two,  and  their  derivatives,  we  can 
use  the  following  properties  given  in  Ref  2,  which  hold  for  integers  v  ^  1  : 


2 


v+1(z) 

m 

—  J  (z)  -  J  ,(z) 

Z  V  v-1 

(15) 

— 

V+1(2) 

a 

f  V2>  -  V,(2) 

(16) 

T 

v+,(z) 

* 

Wz)  -f  Vz> 

(17) 

TT 

:  f  (z) 

V+l 

a 

Kv-1(z)  +  fKv(z) 

(18) 

j;(z) 

m 

i[Jv-t(z)  '  Jv+1(z)] 

(19) 

K(s> 

- 

i[vi(z>  -  Wz>] 

(20) 

Also  the  Bessel  functions  are  related  to  the  modified  Bessel  functions  by  way  of: 


I  (-it) 

v 


(-i)V  J  (z) 

V 


(21) 


K  (-iz) 

v 


i  + 1  Vz)] 


(22) 
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Equations  (3),  (8),  (9),  (12),  (13)  and  (14)  show  that  we  require  K^(tip) 
and  therefore  F±(±ip)  ,  with  p  real  and  positive.  In  equation  (5),  taking  the 
plus  sign  and  putting  s  »  ip  ,  gives  the  singularity  at  w  »  s  in  the  upper 
half  u-plane  so  that  C+  (which  oust  lie  below  it)  can  be  taken  to  be  the  real 
axis;  taking  che  minus  sign  and  putting  s  -  -ip  ,  gives  the  singularity  in  the 
lower  half-plane,  so  that  C_  (which  must  lie  above  it)  can  again  be  taken  to  be 
the  real  axis.  If  p  «  0  ,  the  contour  must  be  indented  at  the  origin. 

From  equation  (4),  K(s)  is  a  real  and  even  function  of  real  s  ,  so 
taking  C±  to  be  the  real  axis  in  equation  (5)  and  taking  upper  or  lower  signs 
throughout,  we  have:  ddd 


F±(±ip) 


(23) 


So  F  +  (ip)  and  F_(-ip)  are  real  and  equal.  We  observe  that  as  w-H)+  , 

Y  “  I <*i|  *  u  .  Then,  using  the  expansions  of  the  Bessel  functions  for  small 
argument  (see  Ref  2)  in  equation  (4) ,  we  find 

iS  K(W)  -  i(l  -  ^2)  • 

Similarly,  using  the  expansions  for  large  argument,  we  find 

llm  wK(u)  -  1  +  \2  ; 

so  the  integral  is  convergent  as  u)-«»  in  equation  (23). 

Equation  (6)  is  now  used  to  evaluate  K  (tip)  .  We  note  that  when  s  ■  ip  , 

l  1 

the  app.opriate  branch  of  (s  +  i6)  (with  positive  real  part)  is: 
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i ir/4.  .  ..J  iu/4  J 

e  (p  +  4)  -*■  e  p* 


as  4-»0  ,  and  when  s  =  -ip  the  appropriate  branch  of  (s  -  i4)!  is 


-iir/4.  ..J  -iir/4  J 

e  (p  +  4)  -►  e  p 


as  4-*0  :  thus 


K±(±ip)  ■  (l  +  X2  j  exp  ^F±(±ip)j  j .  (24 

The  function  iK(p)  ,  given  by  equation  (10),  must  also  be  written  as  a 
real  function  of  the  real  variable  p  for  computing  purposes.  Differentiating 
equations  (21)  and  (22),  putting  v  -  2  and  z  »  p  ,  we  have 

I^(-ip)  -  -  i.J^(p)  (25 


K^(-ip)  -  f|j^(p)  +  iY2(p)]  ’ 


Substituting  equations  (21),  (22),  (25)  and  (26)  into  (11)  and  then  into  (10) 
gives : 

-  j[(<2  -  ')v;  -  i(vi  ‘ <27) 

■  T - -“'fw-1'*)  . 

[(x2  -  ,)v;]  .  [V;  - .  %] 

For  brevity,  we  have  dropped  the  argument  p  of  the  Bessel  functions  J^(p),  etc 
on  the  right-hand  sides.  It  also  follows  from  equation  (25)  that  in  the  integrand 
in  equation  (13)  the  factor  il^dp)  can  be  replaced  by  [-J^Cp)]  • 

2 

'.'.iking  use  of  the  Wronskian  relation 


J2y2  *  J2y2  *  7?  * 


■> :  m,-  ■■ 
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equations  (27)  and  (28)  are  written: 


D( 


-ip)  -  f  (a2  -  i)j2(j*  *  V 


iK(p) 


_ -  2>2-(^)j; _ 

(>2  -  ■)• vi] 2  *  - ')%] 


When  X  ■  1  ,  they  simplify  to: 


D(-ip) 


iK(p)  »  -  irp.J’(p) 


We  remark  that  on  the  regional  interface  r  -  1  ,  x  >  0  ,  the 
boundary  condition,  see  equation  (3-28)  of  Ref  1,  requires 


PH,.,  ■  ‘PHr-,  ' 


This  relation  can  be  verified  from  equations  (12)  and  (13);  we  find 
that 


3  K2(-ipr) 
Jr  D(-ip) 


-  p.iK(p)  , 


and  then  both  sides  of  equation  (34)  are  equal  to  -2\x/(1  +  X^)  + 


2  W  2?  ' 

it/2^1  +  X  )  0 


00 

l  .iK(p).J^(p).K_(-ip).e'PX  .(1 


(30) 

(31) 

(32) 

(33) 

kinematic 

(34) 

on  the  way 

•2  ,  where 

-  Ap)  .  (35) 
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COMPUTATIONS 


4.1  Zeros  of  J' 

For  equation  (3)  we  need  to  calculate  and  store  the  positive  zeros  0  of 

n 

>  which  is  given  by  equation  (19)  with  v  =  2  .  We  find  the  zeros  by 
successive  iteration  using  Newton's  method:  if  oj  is  a  current  approximation 
to  a  zero,  the  next  approximation  is 


Using  equation  (19)  again,  we  get: 

JjCiu)  =  Cto)  -  2J^  (is)  +  Jq(u) 

The  Bessel  functions  J2  ,  are  calculated  from  equation  (15)  with 

v  *  1,2,3  in  succession.  The  iteration  is  taken  to  have  converged  when 
successive  iterates  differ  by  less  than  2.5  *  10  ^  .  At  convergence,  is 

conveniently  available  for  use  in  equation  (3) . 

The  starting  guess  at  the  first  zero  is  w  =  3  .  When  a  particular  zero 

SI  has  been  found,  mindful  of  the  asymptotic  form 
n 

J2(u)  ~  -(^)  sin(u  '  *  'I)  ’ 

we  take  <j  =  0  +  rr  as  the  starting  guess  for  the  next  zei o  SI  ,  .  The 

n  0  s  n+1 

asymptotic  form  of  follows  from  equation  (19)  and  that  of  given  in 

Ref  2. 

To  calculate  typical  values  of  given  by  equation  (3)  on  x  =  0  to 

graphical  accuracy,  we  needed  to  compute  as  many  as  400  zeros  of  for  use 

in  the  series. 

4.2  Evaluation  of  K+(ip) 

For  use  (along  with  equation  (24))  in  the  calculation  of  K+(ip)  for 
equation  (3)  we  have  to  evaluate  the  integral  (23)  for  F+(ip)  at  the  zeros 


12 


p  ■  fl  .  For  use  in  equations  (8),  (9),  (12),  (13)  and  (14),  noting  that 
K  (-ip)  “  K+(ip)  ,  we  must  evaluate  equation  (23)  over  a  range  of  positive  values 
of  p  .  Leaving  the  origin  p  -  0  aside  for  the  moment,  we  assume  p  >  0  ,  we 
break  the  range  of  integration  of  m  into  two  parts  [0,p]  ,  [p,«®]  and 

subtract  the  logarithmic  singularity  at  u  -  0  thus: 


where 


F+(ip)  -  F,  +  F2  +  F3 


P  f  dw  L  wK(id)  ,  oil 

7  /  ln  7772  ln  2 

Jffl+pL  i  +  x  J 


p 

P  f  _ dii 

«  j  2 

i  u>  + 


doj  a) 

2  In  2 

P 


F  -  £  [  In 

J  U)fc  +  p  1  + 


F^,  «  which  now  contains  the  singularity  at  w  *  0  ,  can  be  evaluated  analytically 
& 

(Ref  3)  and  the  result  is: 


F2  '  I1 


1  '“(f)  -  f  ' 


where  G  is  Catalan's  constant,  G  *  0.9159655942  ...  .  Making  the  changes  of 
variable  u  »  pW  for  Fj  and  w  =  p/W  for  F ^  ,  we  get: 


:  »  i  (  dW  ln  2KM1 

1  *  J  w2  +  1  1  * 


I  f  dW 

7  j  w2  * 


_  ln  (p/W)K(p/W) 

1  1  +  x2 
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The  notes  following  equation  (23)  imply  that  the  logarithms  in  Fj  and 
both  go  to  zero  as  W-K)  .  F^  and  F^  are  evaluated  independently,  in  adaptive 
fashion,  by  the  trapezium  rule.  The  values  obtained  by  using  two  and  four 
intervals  in  [0,l]  are  compared,  and  if  they  do  not  agree  to  within  the 
specified  tolerance  TOL  -  10  ^  ,  the  number  of  intervals  is  doubled  until  two 
successive  values  do  agree. 

As  one  might  expect,  it  is  found  that  when  p  is  small,  F^  needs  only  a 
few  subintervals,  and  when  p  is  large,  F^  needs  only  a  few  subintervals  so 
that  no  more  work  is  done  than  is  necessary,  and  also  the  program  does  not  call 
for  values  of  m  so  large  that  the  evaluation  of  K(ui)  fails  when  machine 
limits  are  encountered. 

The  above  procedure  fails  for  p  »  0  .  However,  equation  (B-21)  of  Ref  1 
quotes  the  behaviour  of  K+(ip)  for  small  p  : 

K+(ip)  ♦  (~~ 

This  trend  is  confirmed  by  inspection  of  the  values  of  K+  output  for  small 
values  of  p  .  There  is  therefore  no  need  to  consider  the  integral  equation  (23) 
at  p  «  0  . 

F+(ip)  »  0(1 /p)  (again  con- 


as  p  «°  (36) 

4.3  Evaluation  of  A  from  equation  (7) 

In  equations  (8)  and  (9)  we  first  note  the  behaviour  of  the  integrands 
near  p  «  0  and  as  p-*»  . 

Near  p  -  0  ,  we  can  introduce  the  series  expansions  of  J.  and  Y,  into 

*  2  *  -  ‘ 
equation  (31)  to  show  that  iK(p)  «  0(p  )  and  Jj(p)  *  0(p)  .  Also  K+(ip)  is 

bounded.  Hence  the  integrands  of  C  and  C ,  are  finite  except  for  the  term 

-1  1  4 

0(p  )  of  which  however  is  integrable  and  the  integral  can  be  evaluated 

analytically. 


For  large  values  of  p  ,  from  equation  (23 
finned  numerically),  and  so  from  equation  (24): 


K  (ip) 


■ 


1 


a 
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l  .iK(p).J^(p).K+(ip)  -  o(p'5)  , 

P 

by  using  equation  (36).  The  integral  equation  (8)  for  is  therefore  con¬ 

vergent  as  p-~»  .  Also,  from  Ref  3,  the  mean  value  of  the  right-hand  side  of 
equation  (38)  over  a  period  it  is: 

2 

I  +  X2 

so  that,  for  large  p  ,  the  mean  value  of  iK(p)  .J^p)  .K+(ip)  over  a  period  of 
length  tt  is: 
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(, .  ,*/,* 

Hence,  for  Large  p  ,  in  the  integral  (9)  over  a  period  of  length  tt  the  terms 
0(p  ^ )  in  the  square  bracket  cancel,  the  subintegral  over  such  a  period  is 
0(p  2)  ,  for  suff  iciently  large  P  the  integral  from  P  to  «  becomes  a  sum 
of  such  subintegrals,  convergent  to  a  value  0(P  ^)  ,  and  so  the  integral  (9) 
for  Cj  is  similarly  convergent  as  p-*»  . 

The  above  argument  establishes  that  the  expressions  (8)  and  (9) 
for  and  are  proper.  We  now  turn  to  their  evaluation,  considering 

first.  The  approach  is  to  evaluate  the  integral 


r 

/ 


d£ 

P 


iK(p).J^(p).K+(ip) 


(39) 


for  a  sequence  of  integer  values  of  P  and  to  extrapolate  for  the  limit  P-**»  . 
The  range  [o,P]  is  first  divided  into  the  unit  intervals  [0,l]  ,  [  1 , 2 ]  ...  , 
[p-1,p]  .  The  integral  over  each  of  these  unit  intervals  is  evaluated  by  the 
trapezium  rule,  using  up  to  LNPU  successive  bisections,  but  stopping  when  two 
successive  estimates  of  the  integral  over  the  unit  interval  agree  to  within  a 
specified  tolerance  TOl.  LNPU  and  TOL  are  adjustable  program  constants.  In  this 
way,  with  TOL  set  to  10  ^ ,  the  estimates  of  the  integral  (39)  obtained  with  LNPU 
equal  to  4  and  5  were  found  to  agree  to  five  significant  figures. 

The  effect  of  the  upper  limit  P  was  studied,  bearing  in  mind  that  the 
truncated  portion  from  P  to  infinity  is  expected  to  be  0(P  ^)  .  For  X  *  1  , 
the  following  values  of  equation  (39)  were  obtained: 

P  200  400  800  1600 

C,  -1.430537  -1.489004  -1.530398  -1.559720 

The  limit  of  as  P-*»  can  be  estimated  by  Richardson's  method.  Suppose  we 

are  given  values  of  a  function  f  =  f(p)  for  p  «  P  ,  2P  and  4P  .  Assuming 
— r  7 

an  error  term  0(p  )  ,  where  r  is  an  unspecified  exponent,  r  is  given  by: 


and  then 


2r  -  1 


For  P  •  200  (the  first  three  values  in  the  Table),  we  obtain  2  *  1.41287  ; 

the  last  three  give  2r  »  1.41200  .  So  we  can  assume  as  expected  that  r  -  1/2 
Then,  extrapolating  from  the  last  two  values  with  2r  ■  1.4142  gives: 


-  1  .63050 


We  now  turn  to  as  given  by  equation  (9).  The  integral  is  only  conditionally 

convergent  as  p-*»  ,  and  if  it  were  evaluated  in  the  same  way  as  C.  ,  we  should 
expect  slow  convergence  0(P  >  as  the  upper  limit  P-»»  , 

From  equations  (37)  and  (38),  for  large  p  the  integrand  in  equation  (9)  is 
asymptotically 

(X2  +  1 V  T  2  X2(sin  2p  -  1)  1 

~jL7T7  ' 

The  second  term  in  the  square  bracket  has  mean  value  Q— 2 / (X2  +  1)]  ,  so  the 
whole  square  bracket  has  zero  mean  and  represents  an  oscillatory  function  with 
period  x  . 

The  second  term  in  the  square  bracket  has  the  following  further  properties: 
it  is  non-positive,  and  its  curve  touches  the  p-axis  at  p  *  tt/4  +  Nir  ,  where 
N  is  an  integer.  The  function 


(sin  2p  -  l) 


also  has  all  three  properties  listed  above,  and  indeed  coincides  with  the  second 
term  when  \  «  1  .  Therefore,  for  large  p  ,  the  oscillatory  behaviour  of  the 
complete  expression  (41)  resembles  that  of: 


(>2  *  <) 
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and  we  expect  to  get  better  convergence  as  the  finite  outer  limit  P-*»  if  we 
subtract  this  from  the  integrand  in  equation  (9).  The  subtrahend  can  be 

3 

integrated  exactly  ,  and  the  expression  for  C ^  becomes: 


W-.H 


iK(p)  .Jl(p)  .K  (ip)  +  - t —  (1  -  sin  2p) 

1  l 


('  *  ‘ 1 


-i 


To  remove  the  singularity  0(p  ’)  at  the  origin,  we  subtract  another  term: 


which  has  the  same  behaviour  as  the  integrand  as  p-*»  ,  and  is  exponentially  small 

3 

as  p-*»  .  Again  the  subtrahend  can  be  integrated  exactly  ,  and  the  expression 
for  C2  becomes: 


C 


2 


3^---— 3^  *  I  dp  iK(p).J^(p).K+(ip) 


2P 


(42) 


The  integrand  now  vanishes  at  the  lower  limit,  and  we  can  again  use  the  trapezium 
rule. 

The  integrand  is  still  oscillatory,  and  so  we  evaluate  as  a  succession 

of  loops  of  alternately  positive  and  negative  area  (above  and  below  the  p-axis) . 
The  contributions  from  the  loops  are  stored  separately  in  a  one-dimensional 
array.  A  loop  is  deemed  complete,  and  a  new  loop  is  started,  when  the  integrand 
changes  sign  between  two  integration  points;  the  contribution  to  the  integral 
from  these  two  points  is  split  appropriately  into  a  part  which  is  accumulated 
into  the  now  completed  loop,  and  another  which  starts  the  new  loop.  We  expect 
that  the  sum  of  an  even  number  of  these  loops  will  show  a  more  consistent  trend 
as  P-*»  than  the  result  of  cutting  off  the  integration  at  fixed  values  of  P  , 
which  would  correspond  to  cutting  off  random  fractions  of  either  a  positive  or 
negative  loop. 
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To  apply  Richardson  extrapolation  to  three  partial  sums  of  this  alternating 
series,  we  use  the  loop  counter  L  as  the  variable.  Taking  the  largest  loop 
number  to  be  four  times  the  smallest,  which  should  itself  be  even,  we  stop  the 
summation  at  the  largest  multiple  of  8  loops  allowed  by  the  outer  limit  P  . 
Setting  P  •  800  ,  and  also  putting  X  -  1  which  corresponds  to  the  special 
case  mentioned  in  the  Introduction,  we  obtain  the  results  shown  in  the  Table: 

L  126  252  504 

C2  6.08846  6.14185  6.16526 

These  give: 

2r  -  2.2806,  r  -  1.189,  C2(»)  -  6.1835  . 

Had  we  stopped  short  of  P  «  400  ,  the  Table  would  have  been: 

L  64  128  256 

C2  6.01206  6.08995  6.14276 

which  would  have  given: 

2r  »  1.4749,  r  -  0.561,  C2(»)  »  6.2540  . 

Clearly  Richardson  extrapolation  is  somewhat  uncertain,  though  r  should 
be  heading  for  the  value  1.5  .  However,  a  run  with  P  -  1600  produces  the 
following  Table: 

L  254  508  1016 

C2  6.14231  6.16529  6.13653 

C2  is  not  monotonic  in  this  range,  so  Richardson's  method  fails  altogether. 

This  may  be  a  warning  that  the  asymptotic  behaviour  has  not  yet  been  reached, 
or  it  may  be  due  to  a  lack  of  accuracy  in  the  evaluation  of  C2  .  In  either 
case,  some  caution  is  indicated.  We  propose  for  the  amount  to  take  the 
extrapolated  value  arising  from  P  »  800  above:  C2  *  6.1835  .  From  equations 
(7)  and  (40),  and  this  estimate  for  C2  ,  we  get  a  tentative  value: 

A  -  -  0.2637  .  (43) 
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4.4  The  solution  inside  the  pipe  and  inside  the  jet 

He  can  now  evaluate  the  sum  (3)  for  1  ,  x  4  0  ,  and  the  integral  (12) 

for  the  region  inside  the  jet  r  1  ,  x  >,  0  .  On  the  line  x  -  0  ,  where  the 
exponential  factors  are  all  equal  to  unity,  the  results  are  likely  to  be  least 
accurate,  but  inspection  of  the  computed  results  for  equation  (12)  suggests  that 
for  the  integration  outer  limit  P  «  800  ,  they  are  accurate  to  at  least  three 
significant  figures  so  that  we  can  expect  them  to  be  adequate  for  graphical 
purposes.  The  sunmation  over  the  first  400  zeros  of  in  equation  (3)  should 

give  at  least  similar  accuracy,  since  “  400w  *  1250  . 

On  this  crucial  interface,  the  two  solutions  should  agree.  Fig  2(a)  shows 
their  behaviour  there  for  X  ■  1  ,  and  we  can  see  that  they  do  agree  qualitat¬ 
ively,  but  not  quantitatively;  the  curve  for  ^(x  *  0-,r)  (inside  the  pipe) 
lies  above  that  for  ^(x  *  0+,r)  (inside  the  jet). 

An  exercise  was  mounted  to  see  whether  the  differences  might  be  largely  due 
to  a  small  error  in  evaluating  A  ,  which  is  a  somewhat  uncertain  quantity  in 
view  of  the  investigation  reported  in  section  4.3.  When  (-A)  is  reduced  by  0.01 
to  the  value  0.2537,  a  qualitatively  similar  picture  is  obtained  (Fig  2b)  in 
which  the  differences  are  reversed  in  sign  and  the  curve  for  x  »  0-  now  lies 
below  that  for  x  »  0+  .  At  r  »  0.8  ,  where  the  difference  between  the  two 
curves  is  greatest,  linear  interpolation  of  the  differences  produces  a  new  value 

A  «  -  0.25944  (44) 

and  the  other  values  r  *  0.1 (0.1) 1.0  yield  consistent  values  of  A  which 
differ  by  no  more  than  0.0003  ,  except  for  a  value  (-0.25846)  at  r  ■  0.1  , 
where  the  difference  between  the  two  curves  is  too  small  for  the  interpolation 
to  be  meaningful. 

With  A  given  by  equation  (44),  the  two  curves  for  x  -  0-  and  x  *  0+ 
agree  as  shown  in  Fig  3.  As  a  further  check,  extra  values  for  r  «  0.91(0.01)0.99 
are  plotted  and  also  show  good  agreement.  The  increasing  gradient  towards  r  ■  1 
is  to  be  expected  since  (Ref  1)  the  velocity  field  exhibits  a  square-root 
singularity  at  the  orifice  lip. 

It  is  disappointing  that  the  value  (44)  of  A  required  to  ensure 
a  close  match  of  the  pipe  and  jet  solutions  differs  by  as  much  as  1 .6Z  from  the 
best  prior  estimate,  equation  (43),  of  A  .  The  difficulty  appears  to  lie  in  the 
calculation  of  and  Cj  .  For  instance,  if  we  reject  the  Richardson 
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extrapolation  and  use  the  values  *  -1.55972  and  “  6.13653  obtained  with 
the  largest  values  of  p  (1600)  and  L  (1016)  used,  we  find  A  «  -0.2542  , 
which  is  as  far  above  equation  (44)  as  equation  (43)  is  below  it.  The  discrepancy 
therefore  indicates  that  the  analytic  expressions  given  in  Ref  1  are  not  con¬ 
venient  for  numerical  evaluation,  but  it  does  not  suggest  that  they  are  in  any 
way  inconsistent. 

To  demonstrate  the  relation  between  the  solutions  inside  the  pipe  and 
inside  the  jet.  Fig  4  shows  a  contour  map  of  in  the  region  -1  <:  x  ^  1  , 

0  r  4  1  ,  with  A  given  by  equation  (44) ;  the  contours  vary  smoothly  across 
the  common  boundary  x  *  0  .  The  contour  values  have  been  multiplied  by  100  to 

»  4  .A 

get  three  significant  figures  in  the  standard  plotting  package  used.  No  contours 
appear  at  the  top  left  of  the  Figure  because  the  disturbance  velocity  is  almost 
zero  there.  The  contours  at  the  bottom  right  have  been  omitted  because  the 
velocity  potential  varies  rapidly  there  and  we  wish  to  concentrate  attention  on 
the  interface  region. 

Data  are  supplied  to  the  plotting  routine  at  the  nodes  of  the  grid  shown; 
linear  interpolation  is  used  along  the  grid  lines,  and  the  contours  are  formed 
from  straieht-line  segments  joining  points  on  the  grid  lines.  Near  the  orifice 
lip  this  is  clearly  inadequate,  so  the  solution  there  has  been  computed  at  more 
points  and  plotted  on  a  scale  ten  times  as  large  in  Fig  5.  Again  the  contours 
cross  the  common  boundary  smoothly.  The  representation  is  still  lacking  in 
detail  near  the  lip,  so  a  further  enlargement  by  a  factor  of  ten  is  shown  in 
Fig  6.  Irregularities  now  appear  near  the  interface  at  x  -  0  .  However,  we 
note  that  the  contours  in  Fig  6  are  approximately  similar  in  shape,  indicating 
that  the  3mallest-scale  flow  features  are  now  being  picked  up.  This  scale  is 
remarkably  small  -  the  width  of  the  region  in  Fig  6  is  only  1%  of  the  pipe  radius. 
As  an  indication  of  the  significance  of  the  different  scales,  we  note  the  change 
in  direction  of  the  equipotentials  along  the  bottom  of  Figs  4  to  6.  At  x  »  1  , 
Fig  4,  the  equipotentials  are  approximately  aligned  with  the  pipe  axis;  at 
x  *  0.1  ,  Fig  5,  they  are  approximately  normal  to  the  pipe  axis,  and  at  x  *  0.01  , 
Fig  6,  they  are  again  more  nearly  aligned  with  the  axis.  Along  the  pipe  wall, 
r  »  1  ,  x  <  0  ,  the  equipotentials  are  normal  to  the  wall,  in  accordance  with 
the  boundary  condition,  wherever  the  scale  allows  their  proper  representation. 

We  note,  from  the  relatively  wide  spacing  of  the  equipotentials  towards 
the  top  of  Fig  4,  that  the  disturbance  caused  by  the  external  flow  to  the  flow 
in  the  pipe  has  decayed  to  a  negligible  level  within  one  pipe  radius  upstream  of 
the  orifice. 
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4.5  The  solution  on  the  jet  boundar 


We  turn  now  to  the  solution  on  the  jet  boundary,  represented  in  the  present 
linearized  treatment  by  the  half-cylinder  r  ■  1  ,  x  ^  0  .  It  was  shown  at  the 
end  of  section  3  that  the  kinematic  boundary  condition  is  satisfied  and  that 


-2*  ->3. 

1  ♦  x2  X 


where  G„ 


*Sl(\  +  X2)  0 


00 

/  £  •«<■»• 


J^(p).K_(-ip).e"px  (1  -  Ap)  .  (35) 


This  expression  is  related  to  the  jet  shape,  given  by  F(r,9,x)  «  0  ,  by  the 
kinematic  boundary  condition  7F.7p  •  0  ,  where 

p  »  x  +  e2  iji^Cr.x)  cos  29  +  > 

by  equation  (2),  with  the  axially  symmetric  perturbation  \J>Q  omitted.  If  F  is 
expanded  in  the  form 

F(r,9,x)  -  r  -  1  -  e2  f(x)  cos  29  +  > 

the  above  kinematic  boundary  condition  gives,  to  leading  order  in  e  , 


df  .  y2 
dx  3r 


Hence  the  slope  of  the  jet  surface  in  a  meridian  plane,  9  »  const.  ,  is  given  by 


3r  2  df 

—  *  e  .  —  .cos  29  * 

3x  dx 


2e2  x  cos  29  E  G2  C°S  29 

1  *x2 


to  leading  order  in  e  . 

The  first  term  on  the  right  of  equation  (45)  is  the  leading  term  in  the 
expansion  of  the  solution  given  in  section  3.1  of  Xef  1  for  the  development  of 
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Che  jet  away  from  the  immediate  neighbourhood  of  the  orifice.  When  x  is  large 
enough  for  this  term  to  dominate  the  second,  equation  (45)  represents  a  flatten¬ 
ing  of  the  jet  cross-section,  increasing  with  distance  from  the  orifice,  produced 
by  a  lateral  expansion,  for  8  near  ±ir/2  ,  and  a  fore-and-aft  contraction  for 
9  near  0  and  ±ir  . 

Considering  now  the  second  term  in  equation  (45),  we  note  first  that  it  is 
positive  for  X  «  1  ,  because  the  integrand  in  equation  (35)  is  non-positive  and 
not  identically  zero  for  p  real  and  positive.  The  product  iK(p)J^(p)  0  and 
only  vanishes  at  the  zeros  of  J2  by  equation  (33).  K_(ip)  >  0  by  equation  (24) 
since  F_(-ip)  is  real  by  equation  (23).  The  exponential  is  positive,  and  (1-Ap) 
is  positive  because  A  is  negative.  The  two  terms  on  the  right  of  equation  (45) 
are  therefore  opposed,  and  we  shall  see  that  the  second  term  is  dominant  near 
the  orifice. 

The  integral  in  equation  (35)  closely  resembles  that  in  equation  (12),  so 
G2  can  be  evaluated  by  a  modified  version  of  the  code  written  to  deal  with 
equation  (12),  provided  that  the  integral  is  convergent.  In  fact,  the  change 
from  J2  in  equation  (12)  to  J2  in  equation  (35)  is  unimportant,  but  the  extra 
factor  p  makes  equation  (35)  divergent  at  infinity  for  x  «  0  .  Numerical 
evaluation  is  therefore  restricted  to  x  >  0  .  For  X  «  1  ,  the  variation  of 
G-  with  x  is  shown  in  Fig  7.  We  note  first  that  it  decays  to  zero  as  x 
increases,  so  that  equation  (45)  is  consistent  with  the  form  found  for  the  jet 

away  from  the  orifice.  The  rapid  decay  is  consistent  with  the  proportionality  to 
-3 

x  predicted  by  equation  (3-39)  of  Ref  1  for  the  asymptotic  behaviour  as  x-«°  . 
The  numerical  solution  provides  the  further  information  that,  by  x  »  1  ,  ie  half 
a  diameter  from  the  orifice,  the  second  term  on  the  right  of  equation  (45)  is 
only  1.35%  of  the  first. 

On  the  other  hand,  as  x  tends  to  zero  G2  appears  to  increase  without 
limit,  as  would  be  expected  from  the  comment  above  about  the  divergence  of  the 
integral  in  equation  (35).  Hence,  near  the  orifice  the  second  term  in  equation 
(45)  dominates  and  3r/3x  therefore  changes  sign  on  r  =  1  as  x  increases 
from  zero.  The  interpretation  of  this  is  discussed  later. 

The  calculated  behaviour  of  G2  as  the  orifice  lip,  r  ■  1  ,  x  »  0  ,  is 
approached  can  be  related  to  the  analysis  of  Ref  1  in  the  following  way.  At  the 
end  of  Appendix  B  of  Ref  1,  the  asymptotic  expression  is  obtained: 
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I 


(46) 


valid  for  r  »  1-  as  x-*0+  .  Hence,  for  \  »  1  ,  with  A  given  by  equation  (44) 

G2  ~  0.1035  x“*  . 

Calculated  values  of  are  used  in  the  Table  below: 


X 

0.1 

0.05 

0.02 

0.01 

G2 

0.2215 

0.3736 

0.6661 

0.9801 

4  r 

x  G 

0.0700 

0.0835 

0.0942 

0 .0980 

It  is  apparent  that  the  value  0.1035  from  equation  (46)  is  a  plausible  extrapol¬ 
ation  of  the  calculated  values  of  x^  . 

4.6  Singular  behaviour  at  lip 

Further  analysis  shows  that  the  solution . obtained  in  Ref  1  is  two- 
dimensional  sufficiently  near  the  lip  of  the  pipe,  as  would  be  expected.  The 
singular  behaviour  is  that  of  the  planar  incompressible  potential  flow  round  a 
semi-infinite  plate,  in  which  the  velocity  potential  is  proportional  to 

2  £  y 

o'  sin(x/2)  ,  where  pe  *  is  the  complex  coordinate  relative  to  the  plate  edge 
and  — it  <  x  <  x  ,  so  that  here  pe1*  «  x  +  i(r  -  1)  .  The  streamlines  of  this 
flow  are  sketched  in  Fig  8a.  The  constant  of  proportionality  for  the  flow 
inside  the  pipe  and  jet  follows  from  equation  (46),  so  that 


On  the  inside  of  the  pipe,  x  “  “if  and 


— — — -r-  sin  ^  .  (47) 

1  ♦  A 

0  «  -x  ,  so 

-2-A-7  (--:)*  .  (48) 

1  + 


This  singular  term  makes  no  contribution  to  the  velocity  potential  on  the  jet 
boundary  x  =  3  »  but  a  glance  at  Fig  6  shows  that  varies  noticeably  along 

this  line,  even  close  to  the  lip.  We  therefore  write,  as  an  approximation  to 


* 
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i|<2(1,x)  for  X  •  1  and  for  x  near  zero,  the  function  i|i2A  ,  where 

10°  ?2A  -  2.0728  -  0.4893  (100  x) 

In  this  expression,  the  constant  and  linear  terms  are  a  fit  to  values  of  if<2 
computed  for  r  *  1  and  0  ^  100  x  ^  1  ,  and  the  coefficient  of  the  fractional 

power  follows  from  equation  (48)  with  X  ■  1  and  A  given  by  equation  (44) . 

Fig  9  shows  a  plot  of  the  asymptotic  expression  (49)  with  computed 
values  superimposed,  for  -1  4:  100  x  4  1  .  It  is  not  surprising  that  the  regular 
part  of  the  solution  is  well  represented  by  the  two  regular  terms  in  equation  (48) 
as  shown  by  the  agreement  for  x  >  0  ,  since  the  range  of  values  of  x  is  so 
small.  The  agreement  for  x  <  0  is  less  good,  but  supports  the  validity  of 
equation  (47)  and  the  determination  of  A  for  X  ■  1  .  Outside  of  the  pipe  and 
jet,  an  expression  like  equation  (46)  but  with  a  factor  X  on  the  right-hand 

side  is  quoted  in  Ref  1  for  1 ,  so  a  similar  treatment  would  be  possible  if  we 

were  to  evaluate  the  expression  (13)  and  (14). 

We  have  not  calculated  the  variation  of  A  with  X  .  However,  a  simple 
argument  shows  that  A  is  never  zero.  By  equation  (7),  if  A  is  to  vanish  C1 
must  vanish;  by  equation  (8)  if  Cj  is  to  vanish  the  integrand  of  the  integral 
on  the  right  must  change  sign;  but  by  equations  (23)  and  (24)  K+(ip)  is  never 
negative,  and  by  equation  (31)  iK(p)J2(p)  is  never  positive,  so  the  integrand 
is  never  positive.  Therefore  A  is  non-zero  and  the  flow  singularity  exists  for 
all  values  of  X  .  Ue  recall  its  dependence  on  9  ,  the  angular  position  round 
the  lip.  Omitting  from  equation  (2),  we  have  the  velocity  potential  inside 

the  pipe  and  jet,  close  to  x  •  0  : 


x  >,  0 


2.0700  (-100  x) 


(49) 


x  4  0 


$  »  x  +  i(i2(r,x)  cos  29  +  •  (50) 

The  singular  behaviour  therefore  vanishes  (at  least  to  this  order)  for 
9  «  t  ir/4  ,  t  3ir/4  ,  but  not  elsewhere.  A  similar  singular  behaviour  appears  at 
the  sharp  trail ing-edge  of  a  flat-plate  aerofoil  at  incidence  in  the  absence  of 
circulation;  however,  for  the  jet  no  eigensolution  corresponding  to  the  circulat¬ 
ory  flow  about  the  aerofoil  seems  to  be  available  to  allow  a  Kutta  condition  to 
be  satisfied  at  the  lip. 


TM  Aero  2137 


TM  Aero  2)37 


25 


It  appears  therefore  that  the  singularity  is  an  essential  feature  of  the 
present  treatment.  Since  this  treatment  is  based  on  a  linearization,  it 
represents  a  local  inconsistency  therein:  disturbances  assumed  to  be  small 
emerge  as  being  unbounded.  Moreover,  the  failure  to  satisfy  a  Kutta  condition 
implies  that  the  solution  is  unacceptable,  at  least  locally,  as  a  limiting  form 
of  a  viscous  flow  in  the  limit  of  infinite  Reynolds  number.  The  view  taken  in 
Ref  1  is  that  the  breakdown  of  the  solution  at  the  lip  is  a  local  breakdown  of 
the  linearization  assumption,  and  that  the  difficulty  can  be  resolved  by  the 
solution  of  a  non-linear  problem  near  the  lip.  It  is  expected  that  a  numerical 
solution  would  be  required  and  it  is  not  clear  that  such  a  solution  would  exist. 
Indeed,  Kuchemann^’^  suggests  that  the  inviscid  representation  of  the  mixing  of 
uniform  streams  at  different  total  pressures  downstream  of  a  splitter  plate 
involves  a  convoluted  vortex  sheet,  as  sketched  on  the  right  in  Fig  10,  in  an 
essentially  time-dependent  process.  If  that  were  so,  it  would  have  important 
consequences  for  understanding  and  calculating  the  initial  formation  of  the  jet 
in  cross-flow. 

The  resolution  of  this  difficulty  requires  further  effort.  For  the  moment, 
we  confine  ourselves  to  an  attempt  to  interpret  the  existing  solution  in  as 
meaningful  a  way  as  possible.  From  equation  (50)  we  have,  for  the  circumferential 
component  of  velocity  inside  the  pipe  and  jet,  near  the  orifice: 

r7  ff  ”  "  f  ^  sin  29  +  °(e3)  1  (51) 


Hence,  particles  lying  in  the  two  planes  defined  by  0  «  0  (and  ±tt  )  and  by 
9  ”  ±  ir/2  remain  in  those  planes  and  it  is  instructive  to  sketch  the  streamlines 
they  follow.  We  concentrate  on  the  neighbourhood  of  the  lip  and  the  case  A  =  1  , 
for  which  i|/j  ,  as  given  by  equation  (49),  is  an  adequate  representation  of 
iJ^O.x)  ,  and  consider  9  *  tr/2  .  Then,  writing  B  =  0.1035  ,  by  equations  (49) 
and  (50) 


3x 


r«1 , 8=x/2 ,x<0 


2 

-  e 


0.4893  +  B(-x) 


! 

» 


(52) 


inside  the  pipe.  This  expression  vanishes  at 


-  B 
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for  e  <<  1  .  Since  the  circumf erential  and  radial  components  vanish  by 
equation  (51)  and  the  boundary  condition,  the  point  x  »  xs  ,  r  »  1- ,  6  «  -it / 2 
is  a  stagnation  point  of  the  flow.  We  note  that  it  lies  extremely  close  to  the 
lip.  Downstream  of  the  lip 


3^ 

3r 


r-1 


e»x/2,x>0 


2 

-  e 


r-1 ,x>0 


(53) 


The  contribution  to  the  derivative  on  the  right  is  plotted  in  Fig  7.  The 

flow  is  therefore  inward  across  the  line  r-1  on  a  scale  larger  than  xg  , 
and  we  are  led  to  sketch  a  streamline  pattern  in  the  vicinity  of  the  lip  as  in 
Fig  8b.  On  a  smaller  scale  the  pattern  of  Fig  8a  appears,  and  on  a  larger  scale 
the  global  flow  in  this,  the  lateral  plane,  is  outwards,  as  discussed  in  Ref  1. 

It  must  be  remembered  that  the  present  treatment  is  based  on  the  assumption  that 
the  two  streams,  of  different  total  pressure,  are  separated  by  the  stream  surface 
through  the  lip,  supposed  not  to  depart  too  far  from  the  surface  r  =  1  ,  x  >  0  . 
The  flow  of  ambient  fluid  round  the  lip,  even  on  a  very  small  scale,  as  drawn 
in  Fig  8b,  represents  a  significant  violation  of  that  assumption. 

A  slightly  more  consistent  pattern  is  achieved  by  modifying  the  behaviour 
at  the  stagnation  point,  as  shown  in  Fig  8c.  If  the  streamline  forming  the 
boundary  between  the  two  streams  leaves  the  pipe  wall  tangentially,  the  ambient 
stream  can  stagnate  at  a  pressure  below  the  stagnation  pressure  of  the  jet. 
Although  we  have  not  examined  the  flow  outside  the  pipe  and  jet  in  the  same 
detail,  it  is  to  be  expected  that  a  stagnation  point  will  occur  on  the  outside 
of  the  pipe  at  6-0  (and  it  ) ,  in  the  plane  of  the  cross-flow,  at  a  distance 
from  the  lip  similar  to  -x  .  In  this  case,  illustrated  in  Fig  8d,  the  jet 
fluid,  at  total  pressure  above  ambient,  flows  outward  round  the  lip  and  the 
ambient  flow  stagnates.  On  the  large  scale,  the  flow  in  this  plane  is  inwards, 
as  discussed  in  Ref  1 . 

The  streamline  pattern  of  Fig  8c  is  similar  to  that  on  the  left  of  Fig  10, 

,  ,  ii  5 , 6  , 

which  is  taken  from  Kuchemann  .  In  a  discussion  of  various  ways  m  which  flows 

involving  edge  singularities  can  be  modified  in  order  to  satisfy  a  Kutta  con¬ 
dition,  he  suggests  that  the  unrealistic  jet  flow  on  the  left  might  be  modified 
to  produce  the  flow  on  the  right.  This  involves  shedding  a  train  of  vortices, 
represented  by  a  rolled-up  vortex  sheet,  from  the  lip,  in  an  essentially  unsteady 
flow.  It  is  not  clear  that  the  downstream  growth  of  the  vortices  shown  in  Fig  10 
is  related  to  the  satisfaction  of  the  Kutta  condition.  As  explained  in  Ref  7, 
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the  growth  conforms  with  the  observed  growth  in  the  overall  thickness  of  the 
shear  layer  and  would  require  a  further  mechanism  to  allow  the  vortices  to  grow 
in  strength  as  they  are  convected  downstream. 

We  now  look  at  the  streamline  pattern  on  the  inside  of  the  pipe  just  inside 
the  lip,  ie  on  the  scale  of  xg  . 

On  the  inside  of  the  pipe,  equation  (50),  with  the  approximation  (49) 

for  ,  provides  the  velocity  components  which  define  the  streamlines.  The 

pattern  which  emerges  is  established  in  the  Appendix  and  illustrated  in  Fig  11, 

2  2 

in  terms  of  a  stretched  coordinate,  £  •  -x/5  ,  where  6  ■  Be  .  Fig  11a  shows 

the  pattern  of  the  stretched  streamlines  in  the  limit  of  i-*0  ,  and  Fig  11c  shows 
it  for  an  unrealistically  large  value  6  *  0.4  .  The  similarity  of  these 
pictures  demonstrates  that  the  pattern  of  stretched  streamlines  is  almost 
independent  of  5  for  the  small  values  of  6  to  which  the  present  treatment 
applies. 

Assuming  the  behaviour  of  the  flow  on  the  outside  of  the  pipe  is 
qualitatively  the  same,  we  can  sketch  the  surface  streamlines  very  near  the 
lip  as  in  Fig  12.  Here  the  pipe  has  first  been  cut  along  the  windward  generator, 

0  =  ±ir  ,  then  unrolled  flat,  and  then  unfolded  about  the  lip,  x  »  0  ,  to  show  the 
internal  and  external  flows  together.  The  streamlines  run  smoothly  round  the  lip, 
the  two  flows  colliding  along  a  curve  composed  of  elements  similar  to  the  stream¬ 
line  NS  in  Fig  11c.  This  curve  has  the  behaviour  of  a  three-dimensional  separation 

Q 

line,  as  specified  by  Lighthill  ,  each  element  of  it  starting  at  a  saddle  point 

2 

of  separation  at  -x  ■  $  ,  9  -  0  ,  r/2  ,  ±tt  ,  and  ending  at  a  nodal  point 

of  separation  at  x  «  0  ,  9  »  ±ir/4  ,  ±3ir/4  .  It  is,  of  course,  a  surface 
streamline  of  an  inviscid  flow,  not  a  limiting  streamline  of  a  viscous  flow,  but 
continuity  still  requires  that  the  three-dimensional  flow  leaves  the  surface  in 
its  neighbourhood.  It  thus  takes  over  the  role  of  the  lip  itself  in  the  absence 
of  a  Kutta  condition.  We  must  remember  that  the  scale  of  this  flow  is  extremely 
small  and  that  it  arises  from  a  solution  which  is  not  consistent  on  this  small 
scale . 

Finally,  to  obtain  some  idea  of  the  flow  in  the  pipe  transverse  to  its  axis 
we  may  consider  the  component  of  the  velocity  in  such  a  transverse  plane  as 
defining  a  'secondary  flow'  and  sketch  the  'streamlines'  of  this  'flow'.  We 
must  remember  that  they  are  not,  of  course,  actual  streamlines,  nor  even  the 
projections  of  actual  streamlines.  The  radial  and  circumferential  velocity 
components  are  given,  to  leading  order,  by  using  equation  (50): 


28 


V 

r 


3r  "  6  3r 

i  iS  .  li. 

r  38  "  r 


cos  26 


i|<2  sin  28 


and  so  the  equation  governing  the  'streamlines'  is 


r  4—  ”  rp—  *  -  2  tan  28 

dr 


V(4) 


The  essential  properties  of  i|i.,(x,r)  ,  for  x  =  -0.!  ,  may  be  deduced  from 
Figs  4  and  5,  viz ,  is  negative,  and 


3'l<2 

—  <  0  for  0  <  r  <  rQ 


0  for  r  *  0  ,  rQ  and  1 


>  0  for  Tq  <  r  <  1  , 


where  r^  »  0.75  .  These  properties,  together  with  a  slightly  more  detailed 
examination  of  the  relative  magnitudes  of  ^  an<*  its  derivative  near  r  »  , 

are  enough  to  establish  the  'streamline'  pattern  shown  in  Fig  13.  This  indicates 
that  the  secondary  flow  takes  fluid  from  the  neighbourhood  of  the  plane  of  fore- 
and-aft  symmetry  to  the  sides  of  the  jet. 

5  CONCLUSIONS 

(1)  Ref  1  presents  analytical  results  for  the  flow  of  an  inviscid  incompressible 
jet  emerging  from  a  circular  pipe  into  a  weak  cross-flow.  In  particular, 
expressions  in  the  form  of  infinite  series  or  integrals  are  provided  for  the 
velocity  potential  in  the  four  regions:  inside  the  pipe,  inside  the  jet  and 
outside  both,  upstream  and  downstream  of  the  orifice  plane.  These  expressions 
are  so  involved  that  we  are  unable  to  verify  analytically  the  continuity  of  the 
pipe  and  jet  solutions  across  the  orifice.  Careful  numerical  evaluations  have 
now  verified  the  continuity  of  these  solutions  for  the  case  in  which  the  external 
velocity  component  parallel  to  the  pipe  is  equal  to  the  pipe  velocity  (1  ■  1). 
Confidence  in  the  solution  obtained  in  Ref  1  is  thereby  increased. 


s 
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(2)  Ref  1  showed  that  disturbances  to  the  flow  in  the  pipe  decay  exponentially 
upstream  of  the  orifice.  The  numerical  solution  shows  that  the  decay  is  virtually 
complete  within  one  pipe  radius. 

(3)  The  influence  of  the  orifice  on  the  jet  deformation  decays  only  algebraic¬ 
ally,  according  to  Ref  1.  However,  its  influence  is  down  to  1 .4%  of  the  total 
deformation  within  one  pipe  radius  downstream  of  the  orifice. 

(4)  The  expressions  given  in  Ref  1  are  inconvenient  for  numerical  evaluation. 

(5)  In  Ref  1,  it  is  shown  that  the  solution  obtained  fails  to  satisfy  a  Kutta 
condition  of  smooth  outflow  at  the  pipe  lip  whenever  the  quantity  A  is  non¬ 
zero.  We  show  that  A  is  never  zero  and  that  for  the  case  X  »  1  referred  to 
above  it  is  large  enough  to  affect  the  solution  significantly. 

(6)  Consequently,  the  flow  near  the  pipe  lip  should  resemble  the  two-dimensional 
potential  flow  round  the  edge  of  a  semi-infinite  flat  plate.  It  emerges  that  this 
flow  behaviour  occurs  on  an  extremely  small  scale,  and  that  it  is  embedded  within 
a  markedly  three-dimensional  flow  pattern. 

(7)  If  there  is  no  steady  solution  satisfying  a  Kutta  condition,  a  significant 
source  of  unsteadiness  in  real  flows  at  high  Reynolds  number  is  apparent.  It  is 
therefore  important  to  attempt  to  resolve  the  breakdown  of  the  linearized  approach 
by  examining  the  problem  near  the  lip  on  a  smaller  scale,  as  suggested  in  Ref  1. 
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Appendix 

STREAMLINES  ON  THE  PIPE  SURFACE  JUST  INSIDE  THE  LIP 

The  velocity  potential  0  inside  the  pipe  close  to  x  *  0  is  given  by 
equation  (SO)  and,  introducing  the  numerical  approximation  equation  (49)  to  0, 
we  obtain 

0  -  x  ♦  c2(o. 02073  -  0.4893  x  -  0 . 2070(-x)  ^cos  26  , 

2 

to  order  t  .  Hence  the  dominant  terms,  for  small  x  and  small  t  , 
axial  and  circumferential  velocity  components  are: 

~  2  -i 

30/3x  «  1  +  Be  (-x)  cos  28 

30/r38  »  30/38  »  -  0.4  Be2  sin 

where  B  i  0.1035  ,  as  in  section  4.6.  The  circumferential  component  is  every¬ 
where  small,  so  the  streamlines  are  nearly  aligned  with  the  pipe  axis  unless 
the  axial  component  is  also  small.  The  axial  component  can  only  be  small  if 
cos  28  ^  0  ,  ie,  ir/4  ^  1 0  [  ^  3rr/4  .  It  is  sufficient  to  consider  the  interval 
ir/4  ^  8  ^  it/ 2  ,  since  the  streamlines  in  the  other  8-intervals  (it/ 2  ,  3ir/4) 
and  (-3x/4  ,  — it/4)  of  interest  follow  from  considerations  of  symmetry. 

For  brevity  we  introduce  a  small  quantity 

6  «  Be2  ;  (A-3) 

and,  to  work  on  an  appropriate  length  scale,  we  introduce 

£  -  -  x/62  ,  (A-4) 

noting  that  £  i.  0  in  the  pipe.  Then  the  velocity  components  are 

30/3x  »  1  +  £_*  cos  28  \ 

>  (A-5) 

30/38  *  -  0.4  6  sin  28  ,  ' 


(A—  1  ) 

in  the 
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while  the  streamlines  satisfy  the  differential  equation 


d6 

33e 


3<fr/3x 


(A-6) 


Hence,  along  a  streamline, 


d6 

d5 


0.4 


5 


5*  d3  sin  26 
+•  cos  26 


(A- 7) 


and  we  can  regard  this  equation  as  defining  'stretched  streamlines'  in  the  (0,5) 
plane. 

In  the  (0,5)  plane  we  can  easily  find  the  stretched  streamlines  for  the 
limiting  case  d  ■  0  ,  even  though  we  cannot  return  to  the  physical  (0,x)  space 
for  6  »  0  .  We  recall  that  0  lies  in  the  interval  [tt/4  ,  ir/2]  ,  so 
sin  20  }  0  and  cos  20  ^  0  ,  and  that  6  0  .  Then  the  numerator  of  the 

right-hand  side  of  equation  (A-7)  tends  to  zero  through  positive  values  as  d-K>  , 
while  the  denominator  is  negative  for  5^  ^  -cos  20  ,  vanishes  for  5^  «  -cos  20 
and  is  positive  for  5^  >  -cos  20  .  The  stretched  streamlines  are  therefore  the 
straight  lines  of  constant  0  ,  on  which  d0/d5  is  zero,  together  with  the  curve 
2 

5  »  cos  20  ,  on  which  the  right-hand  side  of  equation  (A-7)  is  indeterminate. 
These  lines  are  drawn  in  Fig  11a,  with  arrows  indicating  the  direction  of  flow 
which  follows  from  equations  (A-5) .  It  is  to  be  expected  that  the  stretched 
streamlines  for  d  non-zero,  but  small  compared  with  unity,  will  resemble  those 
of  Fig  1  la. 

We  now  seek  more  definite  information  about  the  streamlines  for  0  <  6  <<  1 
It  is  helpful  to  locate  and  identify  singular  points  of  streamline  patterns, 
indicated  by  indeterminacy  of  the  expression  for  the  streamline  slope,  and 
usually  associated  with  stagnation  of  the  flow.  Inspection  of  equation  (A-7) 
shows  the  right-hand  side  becomes  indeterminate  for  6  =  ir/2  ,  5  “  1  and  for 
3  *  tt/4  ,  5  -  0  .  Further  consideration  shows  there  are  no  other  singular  points 
Reference  to  equation  (A-5)  shows  that  the  first  of  these  is  a  stagnation  point, 
with  both  velocity  components  in  the  surface  of  the  pipe  vanishing.  Moreover, 
if  we  expand  about  this  point,  writing  0  *  tt/2  -  0  ,  5  «  1  -  5  >  with  0,5 
small  and  of  the  same  order,  we  find  that  d0/d5  ~  -  1.6  d3  0/5  • 


TM  Aero  2137 


Appendix 


33 


From  first  principles,  or  by  consulting  standard  texts,  it  follows  that  the 
point  6  -  x/2  ,  £  -  1  is  a  saddle  point,  at  which  the  streamline  directions  are 
along  8  •  0  ,  8  «  ir/2  and  along  £  ■  0  ,  £  “  1  .  It  is  labelled  S  in  Fig  11c. 
The  second  singular  point  at  8  »  ir/4  ,  £  «  0  ,  labelled  N  in  Fig  11c,  is  not 
a  stagnation  point,  since  3$/38  is  negative  there.  Expansion  about  the  point 
suggests  that  the  local  behaviour  is  essentially  non-linear.  With  8  *  ir/4  +  8  , 
where  8  ,  £  «  1  »  equation  (A-7)  gives 


m  3  i 

d8  0.4  SJ  V 

-  -w  ■>  ■■ 

dS  £  -  25 


(A-8) 


From  this  equation  we  find  that  there  are  again  just  two  directions  along  which 

-2 

strerched  streamlines  reach  the  point  N  ,  namely,  along  £  ■  0  ,  when  £  ~  48  , 

■  3  m2 

and  along  8  «  0.4  6  £  ,  when  £  >>  48  .  The  point  resembles  a  node,  but 

equation  (A-8)  does  not  seem  to  be  equivalent  to  a  linear  system,  so  the  standard 
classification  of  singular  points  may  not  be  applicable. 

We  now  find  a  particular  solution  of  equation  (A-7),  valid  for  5  <<  1  , 

which  connects  the  two  singular  points  S  and  N  .  We  are  guided  by  the  solution 
2 

£  «  cos  28  for  d  -  0  shown  in  Fig  11a  and  repeated  as  the  broken  line  in 
Fig  11c.  For  convenience,  let 

cos  28  ■  -  n  ,  n  >f  0  , 


so  that  equation  (A-7)  becomes 


d£  _  2(1  -  n  £~*  ) 

dn  2\ 

a(l  -  n  ) 


where  a  is  another  small  parameter  defined  by 


(A-9) 


a  =  1.6  53  =  1.6  B3  e6  ,  B  *  0.1035  .  (A-10) 

Then,  if  we  assume  an  expression  of  the  form 


£ 


nZ)(f0(n) 


+  anf 1 (n) 


2  2t/, 
+  a  n  f^(n) 


(A-11) 
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we  have 


n£  *  ”  1  -  $an(i  -  n2)^fg  +  an£l  +  a2r|2f2  +  •••) 

♦3/8t.V|l  -  n2j  ^fQ  +  anf1  +  ...y  -5/l6a3n3(l  -  n2  j  (f Q  +  ...j3 


and 


d?/dn  -  2n  +  an2^3  -  5n2)(fQ  ♦  anf1  +  ...) 


+  an3(i  -  n^fg  +  <*nf|  +  •  ••)  +  an3(i  -  n2j(af1  +  ...  )  . 

Hence,  on  introducing  these  expressions  into  equation  (A-9)  and  equating  the 
coefficients  of  like  powers  of  a  on  the  two  sides  of  the  equation,  we  find: 

a0:  fQ  -  2 

nV..  (3  -  5n2)f0  ♦  n(t  -  n2)f’  -  f,  -  3/4 (t  -  n2)f2 

n3a2:  (3  -  5n2)f,  ♦  n(i  -  n2)f ;  ♦  (t  -  n2)f,  »  f 2  -  3/2(i  -  n2)f0«, 

+  5'8(<  -  St  • 

From  the  structure  of  this  sequence  of  equations  it  appears  that  the  functions 
£q  1  1  fj  1  ...  can  be  found  sequentially  and  explicitly,  and  that  each  is  a 

polynomial  in  r\3  .  In  particular 


f,  -  9  -  I3n2 

f2  -  2(29  -  94n2  +  69n4)  . 
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It  is  not  obvious  that  equation  (A-11)  is  a  convergent  series  in  a  for 
all  n  •  We  therefore  regard  it  as  an  asymptotic  expression  and  re-write  it  in 
the  form 

€  -  ri2  +  aF^n)  +  a2F2(n)  +  a3F3(n)  +  ...  (A-1 2) 

where  F1  ■  2q3(l  -  q2j 

f2  ■  n4(l  -  'i2)(9  -  I3n2) 

f3  -  2n5{l  -  n2)(29  -  94n2  +  69n4)  . 

The  behaviour  of  these  functions  is  illustrated  in  Fig  11b.  Remembering  that 
the  entire  treatment  depends  on  e  being  small  compared  with  unity,  we  regard 
e  »  0.5  as  the  largest  eligible  value.  Then,  by  equation  (A-10),  j  ?  3  *  10  5 
is  the  largest  eligible  value  of  a  ,  so  the  first  two  terms  in  equation  (A-1 2) 
are  sufficient  for  present  purposes. 

In  a  graphical  representation  of  the  stretched  streamlines,  this  particular 

.  .  2 
solution  will  be  indistinguishable  from  the  curve  5  «  cos  20  if  a  is  as 

small  as  3  x  10  3 .  In  order  to  illustrate  the  structure  of  the  stretched  stream¬ 
lines,  an  unrealistically  large  value  of  0.1  is  assigned  to  a  ,  leading  to  the 
particular  streamline  shown  by  the  full  line  connecting  S  and  N  in  Fig  11c. 

We  can  then  sketch  the  stretched  streamlines  between  this  particular  streamline 
and  the  lip  as  shown  in  the  Figure,  noting  that  they  leave  the  lip  with  infinite 
slope,  that  3<j>/ 30  is  negative,  that  3ij>/3x  and  therefore  3$/35  vanish  on 
the  broken  line,  and  that  they  cannot  cross  the  particular  streamline.  The 
stretched  streamlines  upstream  of  the  particular  streamline  can  be  added,  on  the 
basis  that  3<(i/3x  is  positive  and  therefore  3$/3f  is  negative,  that  3$/30 
is  negative  and  much  smaller  in  absolute  value  than  3<f>/ 3 4  except  close  to  the 
particular  streamline,  and  that  they  cannot  cross  the  particular  streamline. 

The  sketch  in  Fig  11c  is  thus  completed. 
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LIST  OF  SYMBOLS 

pipe  radius 

C,/C2  (A  0.26  for  A  -  1) 

constant  (*  0.10  for  A  -  1) 

constants,  see  equations  (8)  and  (9) 

see  equation  (11) 

see  equation  (5) 

see  equation  (35) 

modified  Bessel  function  of  first  kind  and  integral  order  v 

Bessel  function  of  first  kind  and  integral  order  u 

modified  Bessel  function  of  second  kind  and  integral  order  v 

see  equation  (4) 

see  equation  (6) 

see  equation  (10) 

real  variable 

p  -  5ir/4 

distance  from  pipe  axis,  referred  to  pipe  radius  (also  exponent 
in  Richardson's  method) 

complex  variable 

distance  downstream  of  orifice  plane 
x/a 

value  of  x  at  stagnation  point  inside  pipe 
Bessel  function  of  second  kind  and  integral  order  v 

3 

small  parameter  (1.6  5  ) 

2 

Be  (also  arbitrary  small  quantity) 
ratio  of  cross-flow  speed  to  jet  speed 

angular  coordinate  about  pipe  axis,  measured  from  leeward 
generator 

ratio  of  co-flow  speed  to  jet  speed 

2 

stretched  axial  coordinate,  -x/6 


Jk 
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4>(x,r,9) 
$(x,r,8) 
<l«0(r,x) 
iT0(r,x) 
4>2Cr  .x) 
Cr  ,x) 


LIST  OF  SYMBOLS  (concluded) 
velocity  potential  outside  pipe  and  jet 
velocity  potential  inside  pipe  and  inside  jet 
axisymmetric  part  of  $ 
axisymmetric  part  of  0 

2 

non-axisymmetric  part  of  $  ,  to  order  e 

~  2 

non-axisymmetric  part  of  $  ,  to  order  e 
numerical  fit  to  tfi.O.x)  ,  see  equation  (49) 
complex  variable 
n  positive  zero  of 
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Fig  1  Definition  sketch  of  configuration 
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Fig  7 


Fig  7  Comparison  of  normal  derivatives  from  solutions  inside  and  outside  jet  on 
the  jet  boundary 


Fig  8  Streamline  patterns:  (a)  semi-infinite  plate,  (b)  local  behaviour  near  lip  for  e  =  in/2 
(c)  e  =  tn/2  ,  with  different  total  pressures  (distorted  scales),  (d)  as  (c)  for  e  =  0, 


Fig  9  Comparison  of  numerical  and  asymptotic  solutions  on  r  =  1,  very  near  lip 
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